Persistence in higher dimensions : a finite size scahng study 
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We show that the persistence probability P(t,L), in a coarsening system of linear size L at a 
time t, has the finite size scaling form P{t, L) ~ L~^^ f{^) where 9 is the persistence exponent and 
z is the coarsening exponent. The scaling function f{x) ~ x~ for a; <g 1 and is constant for large x. 
The scaling form implies a fractal distribution of persistent sites with power-law spatial correlations. 
We study the scaling numerically for Glauber-Ising model at dimension d = 1 to 4 and extend the 
study to the diffusion problem. Our finite size scaling ansatz is satisfied in all these cases providing 
a good estimate of the exponent 6. 
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Persistence decay has been the subject of considerable 
research activity in recent years. The basic quantity 
under investigation is the persistence probability P{t), 
which is the probability that a given stochastic variable 
with zero mean retains its sign throughout the time in- 
terval [0 : t]. For a large number of systems, it was 
found that at asymptotic times t, P{t) ~ t~^, where 9 is 
in general, a dimension dependent, non-trivial exponent, 
believed to be unrelated to the other known exponents 
W. The non-triviality of 9 is particularly true for spa- 
tially extended systems where the time evolution of the 
stochastic field at one lattice site is coupled to that of 
its neighbours, making the effective single-site evolution 
non-Mar kovian . 

In recent times, the spatial aspects of the persistence 
problem has also come under study. In particular, the 
spatio-temporal evolution of the set of persistent sites 
has been studied by several authors. These include the 
diffusion problem in d = 1 ^J, Ising models in spatial 
dimension d — I |^|-|| and d = 2 Q and the generalised 
g-state Potts model in d = 1 [0. It was found that the 
interplay between persistence decay and the underlying 
coarsening process leads to dynamical scaling and frac- 
tal formation in the spatial distribution of the persistent 
sites. Such fractal structure has also been reported in an 
experimental study of breath figures ||] . 

In the present paper we propose a scaling form for the 
persistence probability P{t, L) as a function of the lattice 
size L and time t. We use physical arguments to moti- 
vate the scaling form in the context of the Ising model 
and show that the scaling reflects the fractal nature and 
power-law correlations in the spatial distribution of per- 
sistent sites. We provide numerical evidence for its valid- 
ity through simulations in spatial dimension d ~ 1 to 4. 
The analysis is further extended to the diffusion problem 
where approximate analytic theories have been used to 
predict 9 in all dimensions. We argue that fractal forma- 
tion in diffusion should take place in all dimensions and 
provide supportive results from simulations. 

Let us consider the Ising model in a d-dimensional ge- 
ometry of linear size L. We start from an initial random 
configuration and quench the system, say, to the tem- 



perature T = 0. As a result, the spins evolve in time 
following the Glauber dynamics, lowering the total en- 
ergy of the configuration in the process. In course of 
time, domains of positive and negative spins form, with 
characteristic length scale £,{t) growing as a power-law in 
time ie., £,(t) ~ t^'^, where z is the dynamical exponent 
for the coarsening process M . The fraction of persistent 
spins decays as power of time : P(i, L) ^ t~^ as long 
as f <C t* ^ i^. For t ^ t* , the domain cannot grow 
any further because of the finite system size and per- 
sistence probability stops decaying, attaining a limiting 
value P{oo, L) ^ L^^-^ . This happens as long as 
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(1) 



For z9 > d, persistence probability will decay to zero 
for any lattice size L. Also we assume that there is no 
'blocking', whereby a finite fraction of spins never flip, 
leading to a limiting value Poo independent of finite size 
effects. Such a situation is believed to occur in Ising 
model for dimensions d > 4 [nO[ and in disordered sys- 
tems [0. 

The above behaviour of the persistent fraction P{t, L) 
for finite lattice sizes can be summarised in the following 
dynamical scaling form. 



P{t,L)=L-^'!{t/L^) 



(2) 



where the scaling function f{x) ~ x~^ for a; <^ 1 and 
f{x) — > constant at large x. Similar finite size scaling 
ideas have been used in a previous work in the context 
of global persistence exponent for nonequilibrium critical 
dynamics [O]. 

The finite-size scaling form given by EqJ2 implies the 
presence of scale-invariant spatial correlations in the sys- 
tem, characteristic of fractals. To show this, we con- 
sider the two-point correlation function C(r,t), which 
we define as the probability of finding a persistent spin 
at a distance r from another persistent spin. For a d- 
dimensional system, C'{r, t) satisfies the normalisation 

condition J^ C{r,t)d'^r — L'^P{t,L). After substituting 
Eq.0, this becomes 



C{r,ty-'dr - L^-^^fit/L') 



(3) 



Let us rewrite this equation in terms of a new function 
F{a, h) = a^^C{a, b) and dimensionless variables x = r/L 
and r — tjU' . 
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Since the RHS of the equation has no exphcit L- 
dependence, LHS should also be likewise. This is pos- 
sible only if F{a, h) = g{ba~^), where 5(77) is given by the 
integral relation 



77^-(i+^).g(ry)d,7 ~ z/(t) 



(5) 



Using the above equation, the limiting behaviour of 
the function g{r]) for small and large values of the argu- 
ment could be deduced from the known behaviour of the 
function fir). Consider r ^ 1, where /(r) is constant. 
From Eq. ^, this implies that g{ri) is constant for large 77. 
In the other extreme of r ^ 1, /(r) ~ t~*. We split the 
integral in Eq. ^ as J^ — /" + J^ and note that g{r]) 
is constant in the second integral for sufficiently large a. 
The second integral vanishes as t~~^ as r — > 0, whereas 
the RHS diverges as t~*. This can be consistent only 
if the first integral diverges as t~^ , which would imply 
that g{ri) ~ 77"* as 77 ^ 0. This leads to the following 
dynamical scaling form for C'{r,t). 



C{r,t) 



-Z0 



9i^) 



(6) 



For small separations r ^ t^/^, this scaling form im- 
plies scale-free correlations, ie., C{r, t) ~ r~^^ , character- 
istic of a fractal with fractal dimension df — d — z9. On 
the other hand, over larger length scales, C{r,t) ^^ t"^, 
which is indicative of the absence of any spatial correla- 
tions. This scaling description was introduced by us pUq] 
in the context oi A + A ^ (l) model, and later verified 
numerically in 2-dimensional Ising model |^ also. 

To check the finite-size scaling form given by Eq. 0, 
we simulate Ising spin systems of various sizes in spatial 
dimension d = 1 to 4. Starting from a random initial 
configuration, the spins are quenched to zero temper- 
ature and are updated sequentially using the Glauber 
updating rule by which a spin is always flipped if the re- 
sulting energy change AE < 0, never flipped if AE > 0, 
and flipped with probabihty i if AE = 0. One MC 
time step was counted after every spin in the lattice was 
updated once. The persistence probability at any time 
t was determined as the fraction of spins that did not 
flip even once till time t since the time evolution started. 
The data is averaged typically over 1000 starting random 
configurations for small L and low d and over 50 starting 
configurations for large L and high d. 




FIG. 1. The persistence probability P{t,L) is plotted 
against time t (measured in MC steps) for three different lat- 
tice sizes L in d = 2 Glauber Ising model. 
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FIG. 2. Same as Fig. 1, except that the scaling func- 
tion f{x) = L^ P{t, L) is plotted against the dimensionless 
scaling variable x — t/L". The data for different L values 
were found to collapse well to a single curve for 9 = 0.21 and 
2 = 2.12 ±0.05. 



For T = Glauber dynamics of Ising model, the persis- 
tence exponent 9 is exactly known to be 3/8 in d = 1 [Q . 
In higher dimensions, simulations predict 9 ~ 0.22{d = 2) 
0,|jjl| and 9 ~ 0.16(d = 3) @. In our finite size 
scaling analysis of the simulation data, we adopt the fol- 
lowing procedure. For d ~ 1,2 and 3, we fix 9 at its 
known value and adjust z to find the value which gives 
the best data collapse. In all cases, we find z ~ 2, which 
is the accepted value of the coarsening exponent for non- 
conserved scalar models [B|. (In d = 3 Galuber dynamics, 
a slower t^'^ coarsening has been observed before p6[ . 
This is presumably due to lattice effects, but we have 
not seen any signature of this effect in our simulations). 
In d = 4, on the other hand, we fix z at 2, and adjust 
9 to collapse the data to a single curve. The results are 
displayed in Figs. 1 to 4. 

In d ^ A, we find that ior z = 2, 9 ^ 0.12 ± 0.02 
gives reasonably good data collapse over the time scales 
and system sizes studied. Fig. 4 shows the scaled data in 
d = 4. It may be mentioned that in d = 4, earlier Simula- 



tions had suggested that the persistence decay might be 
slower than a power-law, and perhaps logarithmic [[lO| . 
However, the agreement of our data with the scaling form 
Eq.g suggests that persistence follows a power-law decay 
in d — A also. For d > 4, blocking of spins has been 
shown to lead to a limiting value of P(i, L) as t ~> cx3, 
which is independent of L ||l^ . Wc could simulate only 
small lattice sizes for d — 5 from which we cannot make 
any conclusive remark at this stage. 
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FIG. 3. The scaling function f{x) = L^^P{t,L) is plot- 
ted against the dimensionless scaled time x = t/L^ for three 
L-values in d = 3 Glauber Ising model. The observed data 
collapse has been obtained for z = 2.05 and Q — 0.166. 
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:r = f/P 
FIG. 4. The figure shows the scaled probability plotted 

against the dimensionless scaled time in d = 4 Glauber Ising 
model. We have fixed 2 = 2, and find that 61 = 0.12 ± 0.02 
gives the best data collapse. 

In the diffusion problem, wc have a scalar field (t>(x, t) 
evolving according to the diffusion equation. The ini- 
tial values (/)(x, 0) are taken to be independent random 
variables with zero mean. 



dt 



(x,t) ; (0(x,O)^(x',O))=<5'^(x-x') 

(7) 



were in good agreement with simulation results. The ex- 
ponent was found to increase with dimension, and has 
been suggested to have the asymptotic value 9{d) ~ a\/d 
as d — i- oo. The constant a has been estimated to be 
~ 0.14 [|l|,|l|l and ~ 0.18 ||l|] by different authors. For 
d — 1,2 and 3, the exponent values are found to be 
e ~ 0.12,0.18 and 0.23 respectively. 

To simulate Eq. m numerically, we use the finite differ- 
ence Euler discretization scheme on cubic lattices of L"^ 
sites 0,|l8|. 



(x, t + At)= (P{yi, t) 



^(/)(x',t)-2d0(x,i) 



(8) 



where x' runs over all the 2d nearest neighbour lat- 
tice sites of X in the cubic lattice and a — t-^At < j^ 
for stability of the discretization scheme. We have taken 
a — 4-; in our simulations as this value has been observed 
to provide the fastest approach to the asymptotic regime 

0- 

For the diffusion problem, simple scaling arguments 
suggest that the dynamical exponent z = 2 in all dimen- 
sions. In all dimensions studied, we found excellent scal- 
ing collapse with z ~ 2 and the 9 values quoted above. 
Upon substitution of the exponent values into Eq. m, it 
can be easily seen that the condition for fractal forma- 
tion is satisfied for d = 1,2 and 3. For d = 1, this has 
already been confirmed by an earlier numerical study [Ej . 
Our results for the persistence probability and the scal- 
ing function for three different lattice sizes in d = 2 is 
displayed in Fig. 5 and 6. 

It is also possible to extrapolate these results to the 
d — !■ oo limit using the asymptotic form suggested for 
9. We see that in this limit, the LHS of Eq. ^ vanishes 
as -y=, leading us to conjecture that fractal formation 
persists in all dimensions for the diffusion problem. 
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For this problem, it has been shown using approximate 
analytic theories [^7|-^that P{t) -- t'" in aU dimen- 
sions. The predicted exponent values in low dimensions 



FIG. 5. The persistence probability P{t, L) is plotted 
against time t (measured as the number of MC steps) for 
three difi'erent lattice sizes L in d = 2 diffusion problem. 
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FIG. 6. Same as Fig. 5, except that the scaling function 
f{x) = L^ P{t, L) is plotted against the dimensionless scal- 
ing variable x — t/L^. The data for different L values were 
found to collapse well to a single curve for 6 = 0.186 and 
2 = 2.05 ±0.04. 

To conclude, we have proposed a finite size scaling 
ansatz for the persistence probability in a coarsening 
system. The scaling form corresponds to the fractal 
structure and dynamic scaling characterising the spatio- 
temporal evolution of the persistent set. We check the 
scaling form numerically for Glauber-Ising model and for 
the diffusion problem. Finite size scaling enables us to 
study persistence reliably in higher dimensions. Our re- 
sults agree with the known values of 6 in the case of Ising 
model(from d = I to 3) and in the diffusion problem (we 
have checked upto d — 3). For d — 4 Ising model, we 
find the signature of algebraic decay of persistence with 
9 ~ 0.12, in contrast with what had been reported earlier 
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